Time-dependent magneto-transport in a driven graphene spin valve 
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Based on the time-dependent nonequilibrium Green's function method we investigate theoretically 
the time and spin-dependent transport through a graphene layer upon the application of a static 
bias voltage to the electrodes and a time-alternating gate voltage to graphene. The electrodes are 
magnetic with an arbitrary mutual orientations of their magnetizations. We find features in the 
current that are governed by an interplay of the strength of the alternating field and the Dirac 
point in graphene: The infiuence of a weak alternating field on the zero bias conductance is strongly 
suppressed by the zero density of state at the Dirac point. In contrast, for a strong amplitude of the 
alternating field the current is dominated by several resonant peaks, in particular a marked peak 
appears at zero bias. This subtle competition results in a transition of the tunnel magnetoresistance 
from a broad peak to a sharp dip at a zero bias voltage applied to the electrodes. The dip amplitude 
can be manipulated by tuning the ac field frequency. 

PACS numbers; 85.75.-d,75.47.-m,81.05.Uw 



I. INTRODUCTION 

The discovery of graphene, a single layer of carbon, 
sparked a burst of research resulting in a variety of fas- 
cinating findingsir— . The modification of the electronic 
properties of graphene under various external conditions 
have also been explored extensively. Particularly inter- 
esting examples are the response to electromagnetic fields 
and the frequency dependence of the conductivity^"—, the 
photon-assisted transports^ and current^, the microwave 
and far-infrared responseiSr— , as well as the plasmon 
spectrumi^"— . In addition, new electromagnetic modes 
akin to graphene have been predicte d^^'^^ and it was ar- 
gued that graphene responds intrinsically in a non-linear 
manner to electromagnetic radiations and may serve thus 
as a natural frequency multiplier with potential applica- 
tions in terahertz electronicslS. 

Recently, considerable attention was devoted to the 
spin-dependent transport in graphene and graphene- 
ferromagnet heterostructures, as prototypical spintronic 
systemsi2r— . Several groups have been successful in 
experimentally contacting graphene to ferromagnetic 
electrodes^S, such as cobalt electrodes^"'^^ and permalloy 
electrodesSS. This rendered possible the demonstration 
of spin injection into a graphene thin film by means of 
nonlocal magnetoresistance measurements. Further ex- 
periments evidenced a rather long spin-flip relaxation 
length fa 1/xTO in a single layer graphene even at room 
temperatures^^. An anomalous cusp-like feature of the 
magnetoresistance versus the applied bias was observed 
in a graphene-based spin-valve devicea^^. Following 
these works. Ding et al^ studied theoretically the spin- 
dependent transport through the graphene spin valve de- 
vice, and pointed out that the cusp-like feature at zero 
bias is due to a subtle interplay of the graphene peculiar 
nature and the conventional spin-valve properties. More 
recently, the effects of the external conditions such as 
disorder— and the magnetic impurities^, on the spin de- 
pendent transport in graphene-ferromagnet heterostruc- 



tures have been addressed. 

On the other hand, the spin-dependent transport un- 
der an alternating (ac) field is relatively less explored. 
In this work, we study theoretically the spin-dependent 
transport through a graphene spin- valve device with stat- 
ically biased ferromagnetic leads having arbitrary spin- 
polarization directions in the presence of an external ac 
field acting on the graphene monolayer in form of a gate 
voltage. The employed method is based on the time de- 
pendent nonequilibrium Green's function approach, as 
described in Refs . 29^30 . It is found that the current 
shows a peculiar behaviour when the strength of the al- 
ternating field is varied. For a small amplitude of the ac 
field, the zero DOS at the Dirac point of graphene sup- 
presses strongly the ac field effect, particularly the peak 
at zero bias is diminished. When the ac field strength 
becomes sufficiently large, the insulator-like properties 
of graphene at the Dirac point are less relevant and a 
prominent peak appears in the differential conductance 
at zero bias. In this case, the TMR exhibits a transi- 
tion from a peak to a sharp dip at a zero bias voltage. 
The dip magnitude can be varied by changing the ac field 
frequency. 



II. THEORETICAL MODEL 

The spin valve device under consideration here consists 
of a single layer of graphene sandwiched between two 
ferromagnetic electrodes. The magnetization of the 
left electrode is assumed to align along the y direction, 
while that of the right electrode deviates from the y 
direction by a relative angle 9. A dc bias V applied to 
the electrodes results in a longitudinal static current in 
graphene. In addition, we consider a harmonic ac gate 
voltage Vg{t) with frequency w and strength Vac being 
applied to the graphene sheet, i.e. 

Vgit) = Vac cos Ljt. (1) 
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The nature of the couphng of Vg (t) to the graphene sheet 
follows from the symmetry group analysis, e.g. as done 
recently in Ref. [3l| in full details. Inspecting the results 
of Ref we conclude that Vg {t) couples to the graphene 
sheet through the Fi representation of the point group 
of D^h at the Dirac points of graphene. In effect, Vg{t) 
acts on the transport as ac Fermi level. Nevertheless, 
as shown below, this simple coupling leads to remark- 
able consequences for the transport and TMR^. In the 
present case the Hamiltonian of this system reads thus 

H = Hl + Hr+Hg + Ht (2) 

where 

Hl=Y^ SkLacl^^CWLa (3) 

Hr = ^{ewR-(rMncos9)clj^^CkRa-Mnsm0cln^Cki^ 

(4) 

where ekao- is the single electron energy associated with 
the momentum k and the spin a in the a — L, R elec- 
trode, cl_^^{ckai7) is the usual creation (annihilation) op- 
erator for an electron with the energy ekao-- 

Having specified the nature of the coupling of Vg to 
graphene we write for the graphene Hamiltonian in the 
tight-binding approximation 

(5) 

where a|g.(ai.(j) creates (annihilates) an electron with 
the spin a on the position of the sublattice A, 
bi cr{f'i,cr) creates(annihilates) an electron with the spin 
a on the position of the sublattice B, tg is the nearest 
neighbor((i, j)) hopping energy in graphene layer. Mea- 
suring the energy with respect to the Dirac point of the 
unperturbed graphene leads to e{t) = Vg. 

In the momentum space Eq.(IS]) can be rewritten as 

q.<r , ^ ^ (6) 
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where (/)(q) = -tg E with Si = f(l,V3,0),J2 = 

1=1 

^{1,-^/3,0), S3 = a(l, 0, 0)(here a is the lattice spac- 
ing). Upon diagonalizing the Hamiltonian (j6]) one finds 
^±(q) = itgl0(q)li which can be linearized around the 
K points of the Brillouin zone with the dispersion given 

by 

E±{ci) = ±VF\q\, (7) 

where vp ~ 3tga/2 is the Fermi velocity of electron. 
From Eq.®, one infers that the ac field couples to the 
two Dirac cones equally as a time-dependent gate volt- 
age. Hence, it does not induce any transitions between 



the valleys. Therefore, one only need to evaluate the 
contributions of the single vallej*^, and then multiply the 
final results by a factor 2. The coupling between the 
electrodes and the graphene is modeled by 

^^^^Y. [^kaqcL.aq. + H.c] , (8) 

kqacr 

where Tkaq is the coupling matrix between the a elec- 
trode and the graphene; N is the number of sites on the 
sublattice A. 

Using the nonequilibrium Green's function method, 
the electrical current can be expressed as 

+C?<(t,ti)]Fi(£)}e— 

(9) 

where Tr is the trace in the spin space, /a(£) is Fermi 
distribution function, Qa{t,t) = G^^g^ ,^,^{t,t') and 

qq' 

g^{t,t') — J2 ^qa q'a(^'*') ^^'^ 2x2 matriccs represent- 
qq' 

ing the retarded green's function and the lesser Green's 
function, respectively. In the calculation of Eq.®, we as- 
sume that the dominant contributions to tunneling stem 
from the electrons near Fermi level, and hence assume 
the linewidth function to be independent of q. Thus, we 
have 

r. ^ ( ri ) (10) 

with r- = 27rE rk*„q7kaq'^(e - eka.). 

k 

To calculate Ga"'{t, t') in Eq.® we carry out the gauge 
transformation 

^^■'^(t,^') =^r(i,t')e"'^''''*'''*'^ (11) 

and substitute in the equation of motion. One obtains 
then 

Q:'''{t,t') = /||g-^.-(e)[l-g^^'^(£)E'-.'^(£)]-ie— 

(12) 

where g-^-{e) = j^Y^fdl^Ae) and Y.^-'^ = TU^Lie) + 
q 

i?^J^(e)i?t] with 5qaV(e) = (e±»r,)^l|0(g)|^ ' ^ud 

^^/cos| -sinf \ 
Vsm^ cos I ; 

Introducing a cutoff kc leads to 

gl{e)^-Fo{e)-z7rpo{e), (13) 

e le2_D2| II 
Fo{e) = ^\n ^ I , p,is) = ^9{D~\s\) (14) 

with D — vpkc denoting a high-energy cutoff of the 
graphene bandwidth, kc is chosen as to guarantee the 
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conservation of the total number of states in the Bril- 
louin zone after the hnearization of the spectrum around 
the K point, this is achieved following the Debye's pre- 
scription. 

The lesser Green function Q^{t,t') can be derived by 
applying the analytic continuation rules (cf. Ref,33) to 
the equation of motion of the time-ordered Green's func- 
tion on a complex contour (Keldysh, Kadanoff-Baym, or 
another choice of contour), 

where g<(ei,£2) = ei(ei)S< (ei, e2)e,"(£2) with Gl'^ie) 
denoting the Fourier transformation of Eq. ([T^ , and 

E<(e,e') = ^2nJ:.U^)U^) 

ran 

x^Lfhi^ + muj)S[€ — e' H- (m — n)uj] 

mn 

xKTrR^ fji{e + muj)6[e — e' -\- {m ~ n)cj] 

(16) 

where Jm is the mth order Bessel function of the first 
kind. The identity e*""'"("*) = Em=-oo •^™(")e™"* is 
used in the calculation of Ea. ([T5)) . 

Substituting Eqs.dll]) and (HSJ in Eq.®, we finally 
obtain the time-averaged current 

I = iEJli^) II^^^{G:{s)RTRR^g:is)Ti^} 

m 

x[/fl(£ + muj) - /l(£ + muj)]. 

(17) 

which is an exact response of the system without im- 
posing any restriction on the amplitude of the exter- 
nal electric field Vac- For a weak ac field {Vac << w), 
Jm(-^) « r(m+i) ('^)"'- ^^^^^ 

tributions to the tunneling from the high sidebands is 
suppressed. While for the large amplitude of the ac field 

(Vac » W), J™{^) ^ ^f^^COsi^ - ^ - f ). Thus, 

more sidebands contribute then to the transport (we note 
however, that we are assuming that the external ac field 
is harmonic. For this reason we cannot describe with the 
present method the case of very short pluses in which case 
the harmonics associated with the pulse width become 
also relevant). In Eq.(fT7|). we further set the symmetri- 
cal voltage divisioni/ZL^ij = Ep ± \eV, and put Ep = 
in the numerical calculations. The TMR can be obtained 
according to the conventional definition 

TMR=M^), ,18) 

where /(O, tt) is the time-averaged current flowing 
through the system in the parallel (antiparallel) configu- 
ration. 



ev^_^=o.6r„ 
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FIG. 1: (color online) The averaged current / as a function of 
the ac frequency for the different ac strength (a), and the ac 
strength for the different frequency (b) for parallel configura- 
tion of the electrodes magnetizations. The other parameters 
are taken as fesT = O.OOiro,P = 0.4, _D = SFo, eV = Fo, 
where Fo stands for the coupling between the scattering re- 
gion and the electrodes. 



III. NUMERICAL ANALYSIS 

Before performing and discussing numerical calcula- 
tions we shall clarify first the assumptions specific to the 
present theory. Adopting the wide bandwidth approxi- 
mation for the graphene spin- valve system we neglect the 
energy dependence of the linewidth functions FJ^(e). Ad- 
ditionally, we assume that the two electrodes are made 
of the same material, thus the degree of the spin polar- 
izations of the left and the right electrodes defined by Pl 
and Ppi can be written as T^j^ — f]^ = Fo(1 ± P) where 
Fq describes the coupling between the graphene and the 
electrodes in absence of an internal magnetization, and 
is taken as the energy scale in the following numerical 
calculations. 

The frequency dependence of the averaged electrical 
current for the different ac field strength is shown in 
FigUJa). Some oscillations of the electrical current with 
the frequency can be seen. These oscillation peaks are 
asymmetric and their magnitudes depend on the weight 
of the different side bands (given by the Bessel functions) 
and is non-universal. When the frequency increases to 
the value hw ~ 0.02eI4c, the frequency dependence of 
the current through the graphene device acts as a pure 
resistance which does not vary with the frequency. Ad- 
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ditionally, one can easily observe that at the low fre- 
quency region, the electrical current increases when the 
ac strength 14c grows. This result is equivalent to the 
one obtained in ferromagnet-insulator-ferromagnet(FM- 
I-FM) system'^'*. The reason for this coincidence is easy 
to understand: the DOS in graphene vanishes at the 
Dirac point. The ac field at the low frequency region 
stimulates the absorption and emission of photons close 
to the Dirac point, and thus can not break the insulator- 
type properties of graphene. In this case, graphene sheet 
can still be viewed as a tunneling barrier similar to FM- 
I-FM system. In the high frequency region, the photon- 
induced channels away from the Dirac point contribute 
predominantly to the transport, thus leading to the in- 
crease of the electrical current with the ac strength. The 
Vac dependence of the electrical current for different fre- 
quencies is shown in Fig. [lib). The current first increases 
with Vac, and then decreases involving small oscillations 
when eVac > 30fiw. The nonmonotonic dependence of 
the electrical current with Vac is different from FM-I-FM 
system. The reason is that the high sideband tunneling is 
dominant in the ferromagnet-graphene-ferromagnct(FM- 
G-FM) system, while the population probability modu- 
lation of the sidebands J^(14c/w) tends to suppress its 
contribution. Hence, their combinations result in the pre- 
diction that there will be a maximum in the current for 
some ac voltage. 
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FIG. 2: (color online) The dc bias dependence of the differen- 
tial conductance G for different ac strength Vac ai fkj — 0.4ro 
(a) and for different frequency to at eVac = Z.TFo (b) for paral- 
lel configuration of the electrodes magnetizations. The other 
parameters are taken as those of Fig[2] 

FigUl shows the dc bias dependence of the differential 
conductance G = dl/dV for the different ac strength 



and frequency ui in the parallel electrodes magnetiza- 
tions. The differential conductance as a function of the 
dc voltage exhibits successive resonant peaks that cor- 
respond to a resonant tunneling through the photon- 
induced sidebands. It is interesting to observe that there 
exists a strong competition between the ac field effect 
and the Dirac point in graphene. For a small ac field 
amplitude, the zero DOS at the Dirac point of graphene 
suppresses strongly the ac field effect, thus diminish- 
ing the central peak (weighted by Jq ) in the differen- 
tial conductance versus the bias, as shown in Figl^l^a). 
When the ac field strength becomes sufficiently large, the 
insulator-like properties of the Dirac point of graphene 
are destroyed, which resembles the applied magnetic field 
case^i. If eV = 0, the ac field still pumps electrons 
through the structure behaving like an effective finite 
DOS. This leads to the appearance of an implicit central 
peak in the differential conductance at zero bias. This 
behavior is a marked difference to the conventional tun- 
neling junction^"—, where the central peak is preserved 
in the entire range of ac field strength. This characteristic 
features suggest that for the graphene tunneling junction, 
it is possible to externally manipulate the central peak 
induced by photon in the conductance by changing the 
ac field amplitude. Additionally, one can find that the 
sizes of the side-band peaks rise monotonously with the 
dc bias voltage for the small ac field amplitude, but de- 
crease however for a sufficiently large ac field amplitude. 
The reason of this behaviour can be traced back to a com- 
bined effect of the ac field and the nature of graphene. 
For the weak ac field, the linear DOS of the graphene 
dominates the transport, and thus modulates the magni- 
tude of each resonant peak. However, when the ac field 
becomes sufficiently strong, the population probability 
of the sidebands is suppressed leading to the decrease of 
the resonant peaks versus the bias. FiglSJ^b) shows that 
with increasing the frequency w of the ac field, the in- 
terval between the resonant peaks increases reflecting an 
increase of the distance between the photonic sidebands. 
The additional phenomenon is that with increasing the 
frequency, the each side-band conductance peak shows a 
slight rise because of the enhance of the side-band con- 
tribution to the tunneling by its population probability 

The dc bias dependence of the differential conductance 
for different temperatures T and angles 9 are shown in 
FigO With increasing the temperature, the resonant 
peaks in the differential conductance decrease and almost 
vanish at larger T. This temperature dependence of the 
conductance peaks is similar to that of a noninteract- 
ing, single-particle resonance in the multi-channel model. 
The mechanism is that at high temperatures not all elec- 
tron states of the low-lying subbands are fully occupied 
due to the occupation of the next subbands'^®. However, 
near V — 0, the conductance is almost independent of 
the temperature. This characteristic feature is different 
from that of the graphene system in the absence of the ac 
field. The zero bias conductance in the latter is sensitive 
to the temperature^. This insensitivity here can be un- 
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FIG. 3: (color online)The bias dependence of the differential 
conductance G for different temperature T at = 0, eVac ~ 
2.7ro,?kj = 0.4ro (a) and for different angle 9 at ksT = 
O.OOiro, eVac = 2.7ro, &j = 0.4ro (b). The other parameters 
are taken the same as those of Fig[2l 



derstood due to a lifting of the insulator-type properties 
of graphene at the Dirac point in the presence of the ac 
field. Fig. ^h) shows that a monotonous suppression 
of the differential conductance with increasing the angle 
6 takes place in the whole dc voltage range. This stems 
from the fact that when 6 changes from to tt the number 
of spin-up and spin-down electrons is rearranged. There- 
fore the couplings for spin-up and spin-down electrons 
become different and the conductance decreases. 

The dc bias dependence of the TMR, defined in 
Eg. ([T8| . for different ac strengths Vac and different fre- 
quencies uj is shown in Fig. |4l At nonzero bias volt- 
age, the TMR has a step-like structure as a function of 
the bias voltage which is caused by photon-assisted ef- 
fects. When increasing the ac field strength, the TMR 
decreases since the subbands induced by the ac field en- 
hance the electron transport through graphene [cf. Fig. 
[Sfa)], however they are spin independent, thus giving the 
same contributions to the spin-dependent transport and 
leading to the decrease of the TMR with the ac strength 
Vac- However for the sufficiently large ac field strength, 
the TMR at the high bias voltage has a slight rise which 
is related to the suppress of the population probabil- 
ity. Remarkably, from FigH^a), one can see that the 
TMR exhibits a salient transition from a broad peak to 
a sharp dip at the zero bias voltage. This is due to the 
strong competition between the ac field effect and the 
Dirac point of graphene. An effective DOS developing 
by the large ac field strength pushes the central region 



FIG. 4: (color online) The dc bias dependence of the TMR for 
different ac strength Vac a.t hoj = 0.2ro (a) and for different 
frequency u at eVac = 3.2ro (b). The other parameters are 
taken the same as those of Fig[2] 

from the insulator-like to a more metal-like behaviour. 
This results in a strong decrease of TMR at zero bias. 
Additionally, one can find that with increasing the fre- 
quency, the steps in the TMR become broad due to the 
increase of the distance between the sidebands, as shown 
in FiglUb). In particular, with increasing the frequency, 
the amplitude of the dip in the TMR at zero bias in- 
creases due to the lift of the zero bias conductance by 
the population probability of the main sideband. 

IV. SUMMARY 

In conclusion, we studied theoretically the spin- 
dependent transport through the FM-G-FM system in 
the presence of an external ac field by means of the time 
dependent nonequilibrium Green's function approach. 
We obtained analytic formulas for the electrical current, 
and found that there exists a strong interplay between 
the ac field effect and the Dirac point in graphene. For a 
small ac field amplitude, the zero DOS at the Dirac point 
of graphene suppresses strongly the ac field effect, and di- 
minishes some of the photon-induced resonant peaks in 
the differential conductance when varied as a function of 
the the bias. For a sufficiently large ac field strength, the 
insulator-like properties of the Dirac point of graphene 
are lifted which leads to the appearance of prominent 
resonant peaks in the differential conductance, particu- 
larly at a zero bias. In this situation, the TMR exhibits 
a transition from a peak to a sharp dip at the zero bias 
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voltage due to this subtle competition mechanism. This 
dip magnitude can be manipulated even by changing the 
ac field frequency. Therefore, it is suggested that for 
a graphene tunnel junction, it is possible to externally 
manipulate the central peak induced by photons in the 
conductance and the zero bias TMR by changing the ac 
field amplitude or frequency. 
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